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^ ■ Abstract. We study the critical frontiers of the Potts model on two-dimensional 

bow-tie lattices with fully inhomogeneous coupling constants. Generally, for the 
^ I Potts critical frontier to be found exactly, the underlying lattice must be a 3-uniform 

• hypergraph. A more general class of lattices are the 4-uniform ones, with unit cells 

contained within four boundary vertices. We demonstrate that in some cases, such 
lattices can be decomposed into triangular cells, and solved using a modification of 
standard techniques. This leads to the exact inhomogeneous Potts critical frontiers on 
' ^ I various lattices, such as the bow-tie lattice with five different couplings, and critical 

Q ■ points for asymmetric bow-tie lattices. 

! 1. Introduction 

in 

00 ■ Since its introduction sixty years ago [H |2], the Potts model has been the subject of 

^ ■ intense study. Beneath its simple definition lies a wealth of interesting mathematical 

problems, many of which remain unsolved. One of these is the determination of the 
^ ■ critical temperature, or more generally the critical frontier, for the transition to the 

disordered phase. Although the problem is solved in one dimension {d = 1), it is 
completely unsolved for d > 3, and in two dimensions, critical temperatures can be 
found exactly only on a certain class of periodic lattices [3]. Here, we aim to extend this 
class. 

The g-state Potts model is defined as follows. Consider a rf- dimensional lattice on 
which we assign to each vertex, i, a spin, cxj, which can be in any of the discrete states 
(1, g). The Hamiltonian for the system is defined as 

-(3H = Kj2Sicr,,,a,) , (1) 

where /3 = l/{kBT), K is the coupling constant and the sum is over all nearest-neighbors, 
{ij). The corresponding partition function is given by 

^ = E^"'"' (2) 
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where {a} indicates that the summation is over all possible spin configurations. 

When A' > 0, it is energetically favorable for neighboring spins to align and we 
are in the ferromagnetic regime of the model. In particular, a ground state (T = or 
K —J- +oo) is one in which all spins are in the same state. Obviously, there are precisely 
q such ground states. When K < 0, we are in the antiferromagnetic regime, and the 
tendency of spins is to anti-align. A ground state {K — )■ —oo) of such a system is now 
one in which neighboring spins avoid alignment as far as the lattice connectivity allows. 
In some cases this competition between anti-alignment and lattice connectivity leads to 
an extensive entropy of the ground state ensemble. 

The formulation of the Potts model given by ([T]) and (I2j) is useful for its connection 
to a physical process. However, it is often more convenient to employ a cluster 
representation. Let us denote the graph (lattice) on which the model is defined by 
G = {V,E), where V is the set of vertices and E = (ij) is the set of edges. Using 
the temperature parameter v = — 1 and summing (|2]) over the spins, the partition 
function on G can be written [1] 



where the sum is now over subsets of edges E', with \E'\ representing the number of 
edges in E' and k{E') the number of connected components (including isolated vertices) 
in the subgraph G' = {V,E'). Note that the anti-ferromagnetic regime corresponding 
to real K is now represented by G [— 1, 0). Even so, the unphysical part of the phase 
diagram, w < — 1, is of some mathematical interest [5] and has been explored by various 
means P [3 Ej for different lattices [HI [HI [ini E] and boundary conditions US] . 

In the ferromagnetic regime {v > 0) the critical frontier is the locus of the 
transition between the high-temperature disordered (paramagnetic) phase and the low- 
temperature ordered (ferromagnetic) phase. Similarly, in the antiferromagnetic regime 
we expect a transition from the paramagnetic phase to a low-temperature phase of 
antiferromagnetic disorder. But in general, in the regime v < the above studies 
[3 El [71 El ISl [ini [m [121 US] reveal a rich structure with several transitions between 
phases with various characteristics. We shall sometimes refer to the combined loci of 
these multiple transitions as the critical manifold. 

From ([3]) , we can see that the limit g — )■ 1 is identical with bond percolation provided 
we set V = p/{l — p), where p is the probability an edge is open in the percolation 
model. This correspondence means not only that any critical point found in the Potts 
model solves the related percolation problem, but that results found in percolation can 
usually be generalized to arbitrary q [3l [T3]. For example, the ingenious argument used 
by Wierman to find the bond percolation threshold of the bow-tie lattice in 1984 [15] 
applies more or less directly to the Potts model, and this was done not long ago by Ding 
et al. [H]. 

In recent work [16], it was shown that the percolation critical surface for the bow- 
tie lattice could be extended to include five different probabilities. This derivation is 
not completely rigorous, even from a physicist's perspective, as it partially depends on 




(3) 
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Figure 1. The bow-tie lattice with five different edge- interactions. 

the use of the negative probabihty parameter regime of the triangular-lattice critical 
manifold. In fact, these parameters could be in the range (—00, 0] and thus were not 
really probabilities at all. Nevertheless, in the final formula all probabilities are positive 
and in the range [0, 1], and there is ample numerical [17] and analytical [18] evidence 
to suggest that the argument is correct. This leads to exact thresholds on many new 
lattices besides the simple bow-tie. 

In this paper, we extend this result to the five-interaction Potts model (Figure 
[1]). The approach, laid out in the next section, is essentially identical to that used 
in [TB]. In the Potts setting, the "negative probability" region corresponds to the anti- 
ferromagnetic regime and thus appears more natural, at least from a physical standpoint. 




(a) (b) 

Figure 2. a) Triangular- type hypergraph on which Potts critical frontiers can be 
found exactly; b) spins on the external vertices. 



2. Exact solutions 

One class of lattices with exactly known Potts critical manifolds is shown in Figure [2^, 
where the shaded triangle stands for a generic configuration of spins and edges. In terms 
of the spins on the external vertices (Figure [2]d) this corresponds to effective interactions 
between groups of two and three spins. These interactions can be expanded in terms 
of clusters, generalizing [19] the construction of ([3]). The critical point of this system is 
then given by [3l [19] 

qA-C = 0, (4) 

where A is the weight of a triangle inside which none of the external spins are connected, 
and C is the weight for a triangle with all three spins connected. Applying this formula 
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to the simple triangular lattice, with the edge interactions shown in Figure ^1^, we find 
A = 1 and C = uvw + uv + uw + vw leading to the critical frontier 

q — uvw — uv — uw — vw = . (5) 

Removing the w-edge by setting w = 0, we get the critical frontier for the square lattice 

q — uv = Q . (6) 

In equation ([5]), we can single out an edge, say w, and write 

q-uv 

w = . (7) 

uv + u + V 

This form makes clear the content of the solution; the couplings u and v can be chosen 
arbitrarily, and then ([7]) provides the critical value of w. Note the presence of the square 
critical frontier, equation ([6]), in the numerator of ([7]). The square lattice is super-critical 
if uv > q, and thus on the triangular lattice if we choose uv > q the lattice is super- 
critical already and the role of the w-bond in equation ([7]) is to return the system to 
the critical point. This is accomplished by assigning w an anti- ferromagnetic coupling, 
i.e., ~l <w <Q. 





(a) (b) 



Figure 3. a) Triangular lattice with three different Potts interactions. The critical 
frontier is given by equation ([5]); b) generalization to six interactions. The critical 
manifold is unknown in general, but imposing the conditions ([5]) and ^ gives a critical 
system. 



A generalization of the inhomogeneous triangular system is that shown in Figure 
|3b, in which we have two types of triangles, (ui, fi, wi) and (u2, t'2, W2)- We can obtain 
a critical system by independently setting 

q — uiViWi — uiVi — uiWi — viWi = (8) 

and 

q - U2V2W2 - U2V2 - U2W2 - V2W2 = . (9) 

That this system is critical is a consequence of the difference property of quantum 
integrable systems (see the Appendix for the technical details). Of course, this is not 
the full critical frontier of the 6-interaction triangular lattice, which is unknown, rather 
these formulas occupy a sub-region of the full manifold. 
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Figure 4. A 3-uniform hypergraph on which Potts critical frontiers can be found 
exactly. Triangles of the same color have the same interactions on the edges. 




Figure 5. a) Bow-tie lattice with three interactions for which the critical frontier is 
identical with that of the simple triangular lattice, equation ([S]); b) generalization to 
six interactions. 



3. Bow-tie lattice 



The critical manifold ([5]) follows from a duality argument, the main ingredient of which is 
the fact that the unit cell of the triangular lattice lies between three boundary vertices 
as does that of its dual, the hexagonal lattice. In addition to lattices in which the 
triangular faces are more complicated [3], there are also alternative arrangements of 
triangles that have this property and thus their critical manifolds can be found exactly. 
An example can be seen in Figure H] (ignoring the colors for now). If the triangles are 
given the interactions {u, v, w) as in Figure E^, the critical frontier for this lattice is given 
by ([5]), i.e., identical to that of the ordinary triangular lattice. Now, the two to-edges in 
parallel may be combined into a single one, z, using the parallel reduction formula [20] , 

z = w^ + 2w (10) 

and we may substitute equation ([7]) for w to get the Potts critical manifold for the 
bow-tie lattice, 

q{q + 2u + 2v) - 2u^v - 2uv^ - u^v^ 
- z{u^ + + 2uv + 2u^v + 2uv'^ + u^v"^) 

= (11) 

generalizing to arbitrary q the well-known result for percolation [15] . 

The formula f lTTj) was derived in a recent work of Ding et al. [H]. However, 
somewhat more is possible. Here, we show that the inhomogeneous interactions can 
be extended over neighboring triangles, giving a critical frontier involving all five 
interactions in a single formula. To do this, we note that, as for the ordinary triangular 
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lattice, it is not necessary that all triangles be identical. The left- and right-pointing 
triangles may have different interactions, as indicated by the colors in Figure with 
the assignments shown in Figure |5b. As before, we can set the blue and red triangles 
separately to their critical values by requiring 



and 



q-viui 

wi = (12) 

UiVi +U1+V1 

q-V2U2 , . 

W2 = ^ ^ , 13 

U2V2 + U2 + V2 

and the whole system is therefore at a critical point. However, we may now use the 
parallel reduction formula [20] for the abutting ty-interactions to combine them into the 
effective z according to 

z = W1W2 + W1+W2 . (14) 

Now we simply substitute (|T2l) and (|T3i) into (IT^ and multiply away denominators to 
get 

q{q + Ui + U2 + vi + V2) - U1U2V1V2 

- U1U2V1 - U1U2V2 - U1V1V2 - U2V1V2 

- z{uiU2 + U2V1 + U1V2 + V1V2 + U1U2V1 + 

U1U2V2 + U1V1V2 + U2V1V2 + U1U2V1V2) = (15) 

which is the critical manifold of the bow-tie system shown in Figure [H This reduces to 
the previously-obtained critical surface for percolation [HI [16] by setting q = 1, but is 
a new result for the Potts model. It reduces to the formula ( ITT]) of Ding et al. in the 
special case Ui = U2 = u and f 1 = f 2 = v, and to Wu's checkerboard formula [21] if 
z = 0, providing another derivation of that result. 

A word is in order about the meaning of the equation f lTSj) . Although we represent 
the effect of the two interactions Wi and W2 by a single interaction, z, both the w-bonds 
are really still in the problem. Consider the checkerboard = 0. Here, we choose 

Ml and Vi arbitrarily, and then wi must take its critical value (fT2|) . Then, we select 
W2 according to iHM in order to enforce z = 0. We can now choose an arbitrary value 
for U2, but V2 must be determined from iQ to enforce criticality of the left triangle in 
Figure |5b. The result is that we have chosen four interactions arbitrarily and were able 
to find the value of the fifth to give a critical system. This procedure is represented by 
the single formula ( !T5|) . From this discussion, we can see that the statement "2; = 0" 
really corresponds to a family of values {wi,W2) that exactly cancel this interaction, 
from which we must draw one member according to the values of the other interactions 
in the problem. This reasoning extends to any other value chosen for z. 



4. Asymmetric bow-tie lattices 



Aside from assigning different probabilities to the right and left sides, we may use more 
complicated triangular faces in place of the shaded red and blue triangles in Figure IH 
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(a) 



(b) 



Figure 6. Using a) for the blue (left-pointing) triangles in Figure |3^, and b) for the 
red triangles leads to the Potts critical frontier of the lattice in Figure [TK- 



and consequently, it is not necessary that the two sides contain the same graph. We 
may, for example, use the graph in Figure for the blue triangles in Figure |4|l, and 
the graph in Figure [Hb, which is just the unit cell of the martini-A lattice for 
the red. The critical frontier for simple triangle is given as before, 

q-xy 

W2 = ■ — . (16) 

xy + X + y 

To find the critical point of Fig. [6]d, we employ the condition ([5]) with A, the weight of 
configurations in which none of the boundary vertices are connected, given by 

A = q'^ + q{r + s + t + u + v) + st + rt + rs 

+ tu + tv + uv + ru + sv + rst (17) 

and C, the weight in which all three are connected, by 

C = rstuvwi + rstuv + rstuwi + rstvwi 
+ rsuvwi + rtuvwi + stuvwi + rsuv 
+ rtvwi + rsvwi + ruvwi + rsuwi 
+ stuwi + ruvwi + rtuv + stuv 

+ stvwi + rtuwi + qirvwi + suwi), (18) 

with the result 



where 



= ^ (19) 



a = ([^ + (j^{r + s + t + u + v) 
+ q{rs + rt + st + ru + tu + sv + tv + uv) 

+ qrst — rsuv — rtuv — stuv — rstuv (20) 



and 



fi = q[su + rv) + rsu + rtu + stu + rsv 
+ rtv + stv + ruv + suv + rstv + rstu 

+ rsuv + rtuv + stuv + rstuv . (21) 

Now we plug ( fT6|) and ( fT9|) into ( fT4|) and multiply by I3{x + y + xy) to get the full 8- 
interaction critical frontier. Setting all interactions equal gives the homogeneous critical 
frontier displayed in Table [IK- 



Potts critical frontiers of bow-tie lattices 



8 



Many constructions of this type are possible and we give a very small sample 
in the rest of Figure U\ with their homogeneous critical frontiers also reported in 
Table [H The inhomogeneous frontiers for all the lattices of Figure [7] are included in the 
supplemental material to this submission in the file SJ12_EPL.m, which can be processed 
in Mathematica or opened with an ordinary text editor. 

The argument that leads to these results applies equally well to the configuration 
in Figure Hb, corresponding to rotating every second row by 180 degrees. Again using 
the graphs of Figures and b for the blue and red triangles gives the lattice in Figure 
[8^. The other lattices in Fig. |8] are similarly obtained from those in [7] and have the 
corresponding critical frontiers. 




(c) (d) 
Figure 7. Asymmetric bow-tie lattices. The critical frontiers are given in Table [TJ 




Figure 8. Asymmetric bow-tie lattices with the same critical frontiers as those in 
Figure [T] 
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lattice 


frontier 


a 


q{q^ + 

— Zq U 


7q\ + 20gw^ + 24w^ + qv^) 

'iriiiS 0^n,5 Oc;n,6 Q„,7 „,8 


b 


q{q^ + 
-29v^ 


8qh! + 26qv^ + 35v^ + 2qv^ + 5w^) 
- 24:V^ - 8v'^ - 


c 


q{q^ + 
-Uv^ 


7q\ + 2Qqv'^ + 2Av'^ + gt;^ + 2^;^) 
- 7v'^' - v'^ 


d 


qiq"^ + 
+64^;^ 
-39v^ 


9q^v + 35g^t;^ + Tlgt"^ + 2q^v'^ 
- 2qv^ - 31t;^ - - 



Table 1. Critical frontiers for the lattices shown in Figures [7] and HI 



5. Discussion 

We have presented new exact critical frontiers for the Potts model on certain two- 
dimensional lattices. These new solutions employ the special structure of the bow-tie 
lattice, enabling connection of neighboring triangles and leading to exact solutions on 
lattices that are realizations of 4-uniform hypergraphs. 

Note that these lattices are not of the typical exactly solvable type in which the 
unit cell is contained within three boundary vertices, or terminals; that is, they are 
not realizations of 3-uniform hypergraphs. Rather, these are /our-terminal lattices. 
However, they are a very special subset of this class in that they can be derived from 
three-terminal lattices by fusing two edges. This fact leads to the related observation 
that the duals of these lattices are also four-terminal lattices, a property not generally 
shared by such graphs. Our method is thus not a recipe for finding the critical frontiers 
of general four-uniform hypergraphs, such as the kagome lattice. Its dual, the diced 
lattice, is not of the four-terminal type. 

We finally mention that we have verified our results by computing the critical 
polynomials [8] corresponding to the lattices of Figure [7l using square bases of size 
2x2, and Figure [H] using rectangular bases of size 1x2 unit cells. We found that 
for each of these lattices, the critical polynomial contains a factor which coincides with 
the frontier given in Table [H This is in line with the observations of [HI HI] that exact 
solvability implies that the smallest possible basis (1x1) produces the exact critical 
frontier, and that larger bases (such as 2 x 2) factorize, shedding a factor which is equal 
to the exact critical frontier. 
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(b) 



(c) 




Figure 9. Potts interactions on the lattices in Figures [7] and [H 




Figure 10. Part of the medial graph M.{G) corresponding to a left-pointing triangle 
of G. 



Appendix: Relation to quantum integrability 



It is well-known [2l] that the Potts model on any planar graph G is equivalent to a 
loop model on the associated medial graph }A[G). The portion of the medial graph 
corresponding to the leftmost triangle in Figure |3^ is shown in Fig. [TOl 
The integrable R-matrix of the loop model takes the form [25] 

^,(^) = sin(7-^)l, + sin(^)E,. (22) 

We interpret the coefficients as the set of Boltzmann weights associated with a vertex i of 
J\A{G) (i.e., with an edge e of G). Here 6 is the difference of the two spectral parameters 
carried by the oriented lines of M.{G) that intersect at i. The Temperley-Lieb generator 
Ej can be interpreted graphically as the contraction of the two lines going into and 
the identity operator Ij is the reflexion of those lines at i. The equivalence between the 
Potts model and the loop model implies that the weight of a loop is n = ^/g, and we 
have parameterised this as n = 2 cos(7) in ( l22l) . Moreover, the edge weights u get traded 
for new weights x = uq~^^'^ that are ratios of those appearing in front of the terms in 
( 122|) . More precisely, we have x = sin(6')/ sin(7 — 6) if the two oriented lines of M{G) 
go into the edge e from the same side, and x = sin(7 ~ 6)/ sin(6') if they go into e from 
opposite sides. 

Returning to Figure [10] we note that integrability requires that the spectral 
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parameters r, s and t follow the lines of Ai{G) through the intersection at vertices. 
In particular, on the right side of the figure the s-line is on the top and the t-line is 
on the bottom. But since in the bow-tie lattice the left-pointing triangle is followed by 
a right-pointing one (see Figure [3]), the two lines get permuted back after traversing a 
complete bow-tie motif. In particular, the assignation of spectral parameters to each 
"horizontal" line is compatible with periodic boundary conditions. 

For a choice (r, s, t) of spectral parameters, the edge weights in the Potts model are 
therefore 

u = q^'^ 



sin(s 


— r) 


sin (7 — 


s + r) 


sin(7 — 


t + r) 


sin(t 


— r) 


sin(t 


-s) 


sin (7 — 


t + s) 



w = q'l' — — ^ ^— . (23) 

sm(7-t + s) ^ ' 

It follows that the quantity uvw + uv + uw + vw — q appearing on the left-hand side of 
([5]) is proportional to 

2 cos(7) sin(s — r) sin(7 — t + r) sin(t — s) 
+ sin(s — r) sin(7 — t + r) sin(7 — t + s) 
+ sin(s — r) sin(t — s) sin(t — r) 
+ sin(7 — t + r) sin(t — s) sin(7 — s + r) 

— sin(7 — s + r) sin(t — r) sin(7 — t + s). (24) 

But this is identically zero, by a trigonometric identity. The condition (|5]) is therefore 
simply equivalent to the integrability of the Potts model. The same reasoning applies 
of course to the right-pointing triangles, after some relabeling of the parameters. In 
particular we infer that transfer matrices T(s, t\r) and T(s', t'\r) with different choices of 
the horizontal spectral parameters commute, and hence the transfer matrix commutes 
with an infinite set of conserved quantities. It follows that the bow-tie lattice Potts 
model satisfying (IH])-(|9]) is critical. 

Summarizing, the appearance of the constraint ([5]) between the edge weights 
{u,v,w) is due to the fact that the three spectral parameters {r,s,t) admit only two 
independent differences r — s and s—t (since r—t = {r — s) + {s — t)), and the requirement 
that the integrable R- matrix R{0) satisfies the difference property (i.e., its argument 9 
is the difference of two spectral parameters). 
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